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Abstract 

Bootstrap percolation has been used effectively to model phenomena as di- 
verse as emergence of magnetism in materials, spread of infection, diffusion of 
software viruses in computer networks, adoption of new technologies, and emer- 
gence of collective action and cultural fads in human societies. It is defined on 
an (arbitrary) network of interacting agents whose state is determined by the state 
of their neighbors according to a threshold rule. In a typical setting, bootstrap 
percolation starts by random and independent "activation" of nodes with a fixed 
probability p, followed by a deterministic process for additional activations based 
on the density of active nodes in each neighborhood (8 activated nodes). Here, 
we study bootstrap percolation on random geometric graphs in the regime when 
the latter are (almost surely) connected. Random geometric graphs provide an 
appropriate model in settings where the neighborhood structure of each node is 
determined by geographical distance, as in wireless ad hoc and sensor networks 
as well as in contagion. We derive bounds on the critical thresholds p' c ,p" such 
that for all p > p"(6) full percolation takes place, whereas for p < p'Jff) it does 
not. We conclude with simulations that compare numerical thresholds with those 
obtained analytically. 



1 Introduction 



Some crystals or lattices studied in physics and chemistry can be modeled as consist- 
ing of atoms occupying sites with specified probabilities. The lattice as a whole would 
then exhibit certain macroscopic properties, such as (ferro)magnetism, only when a 
sufficient number of neighboring sites of each atom are also similarly occupied. In 
computer memory arrays each functional memory unit can be considered as an occu- 
pied site, and a minimum percentage of functioning units are needed in the vicinity of 
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each memory unit in order to maintain the array with proper functioning. In adoption 
of new technology or emergence of cultural fads, an individual is positively influenced 
when a sufficient number of its close friends have also done so. 

All three examples cited above may be modeled via a formal process called "boot- 
strap percolation" which is a dynamic process that evolves similar to a cellular automa- 
ton. Unlike cellular automata, however, this process can be defined on arbitrary graphs 
and starts with random initial conditions. Nodes are either active or inactive. Once 
activated, a node remains active forever. Each node is initially active with a (given) 
probability p. Subsequently and at each discrete time step, a node becomes active if 9 

of its nearest neighbors are active, for a fixed value of 9 = 1, 2, 3, As time evolves, 

a fraction of all the nodes are activated. The emergence of macroscopic properties 
of interest typically involve <3? to be at or close to 1 . 

Gersho and Mitra [12] studied a similar model for adoption of new communica- 
tion services using a random regular graph and obtained (implicit) critical thresholds 
for widespread adoption. Chalupa et al |@] were the first to introduce bootstrap per- 
colation formally to explain ferromagnetism. Their analysis is carried out on regular 
trees (Bethe lattices) and a fundamental recursion is derived for computation of the 
critical threshold that has since been used extensively. In the more recent past, results 
for non-regular (infinite) trees have also been derived by Balogh et al |5[]. Aizenman 
and Lebowitz [1] studied metastability of bootstrap percolation on the ci-dimensional 
Euclidean lattice Z d which has now been thoroughly investigated in two and three di- 
mensions, see talnD. The existence of a sharp metastability threshold for bootstrap 
percolation in two-dimensional lattices was proved by Holroyd and recently gen- 
eralized to (i-dimensional lattices by Balogh et al [4]. Even more recently, bootstrap 
percolation has been studied on random graphs G(n,p) by Luczak et al jl9n . In j25ll 
Watts proposed a model of formation of opinions in social networks in which the per- 
colation threshold is a certain fraction of the size of each neighborhood rather than 
a fixed value, a departure from the standard model that is used by Amini in (2J] for 
random graphs with a given degree sequence. 

Many diffusion processes of interest have a physical contact element. A link in 
an ad hoc wireless network, a sensor network, or an epidemiological graph connotes 
physical proximity within a certain locality. Study of diffusion of virus spread in ad 
hoc wireless, sensor or epidemiological graphs requires this notion of neighborhood 
for accurate estimation of likelihood of full percolation. This is in contrast to models 
with long-range reach where physical proximity plays little, if any, role. The natural 
random model for such phenomena is the random geometric graph. In this work, we 
focus on bootstrap percolation on random geometric graphs, a topic that has not been 
investigated, to the best of our knowledge, and obtain tight bounds on their critical 
thresholds for full percolation. 
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2 Random Geometric Graph Model 

One of the transitions from the random graph model G(n,p) of Erdos and Renyi ifToL 
and Gilbert lfl3ll to models that may describe processes constrained by geomet- 
ric distances among the nodes is the model of random geometric graphs (RGGs) by 
Gilbert 11411 . The RGG model has been used in many disciplines: for modeling of 
wireless sensor networks |23|], cluster analysis, statistical physics, hypothesis testing, 



spread of computer viruses in wired networks, processes involving physical contact 



among individuals, as well as other related disciplines, see M22H for more details. For 
example, a wireless sensor network typically contains a large number of randomly de- 
ployed nodes with links determined by geometric proximity enabled by (a small) radio 
range among the nodes that is sufficient to enable successful signal transmission across 
the network. A further application of RGG is in representing d-attribute data where 
numerical attributes are used as coordinates in R. d and two nodes are considered con- 
nected if they are within a threshold (Euclidean) distance r of each other. The metric 
distance imposed on such a RGG captures the similarity between data elements. 

Consider an RGG in two dimensions that is constructed by drawing n nodes uni- 
formly at random within [0, l] 2 and connecting every pair of nodes at Euclidean dis- 
tance at most r. Let us denote this process by RGG(n, r). A summary of basic struc- 
tural properties of RGG(n, r) is as follows. 

(i) RGG(n, r) is a 'homogeneous' geometrical model where the distribution of the 
number of nodes within a distance r from a given node follows the same bino- 
mial distribution Bin(n — 1, r 2 ir) (with appropriate correction when the center is 
within a distance r of the boundary). The average degree D = E(deg) of a node 
is nr 2 ir in the limit. 



(ii) There is a critical value A c such that for r > a/ X c /n there exists a giant com- 
ponent, i.e., the largest connected component of order 0(n) nodes contained in 
RGG(n, r) whpfl [22]. We denote the critical threshold for existence of a giant 
component by r c := y 7 A c /n. 

(iii) In this regime, the second largest component is of order C(ln 2 n). 

(iv) The exact theoretical value of the constant A c is not known. It is experimen- 
tally established that A c « 1.44 for the dimension d = 2 [24], while theoretical 
bounds A c G [0.696, 3.372] are given in [20]. There has been a recent improve- 



ment of the lower bound A c > 4/(3>/3) « 0.7698 



'Whp or "with high probability", means with probability one as n, the number of nodes, tends to 
infinity. 
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(v) RGG(n, r) is connected whp for r > ^Inn/im, 1 1(1 21]. We denote the critical 
threshold for connectedness by r t := yj\a.n/im. 

(vi) Every monotone property in a RGG(n, r) (e.g., existence of a giant component 
and connectedness) exhibits a sharp threshold 11511 . 



In order to simplify our analysis on RGGs, we now introduce G n ^ r which is asymp- 
totically isomorphic to RGG(n, rn -1 / 2 ). Let X be a Poisson point process of intensity 
1 on K 2 . Consider points of X contained in [0, y/n\ 2 representing the nodes of a graph 
denoted G rhr . Two nodes of G njT are connected if their Euclidean distance is at most 
r. Our analysis from here on will be based upon the fact that an instance of G n r is 
isomorphic to an instance of RGG(ra, rn" 1 / 2 ) whp 12211 . 



We parameterize r = W _1 alnn by introducing a new parameter a which mea- 
sures how denser G n , r is compared to an instance G n n at the threshold for connected- 
ness rt. The condition a > 1 enables us to deal with an asymptotically connected 
G n ,r 11161 l22ll . Notice that for sufficiently large n the expected degree is concen- 
trated around its mean a In n, which can be easily derived from the Chernoff and union 
bounds. 

For n = 1000, the critical thresholds for the existence of a giant component and 
connectedness in G n>r satisfy r c w 0.0316 and r t ~ 0.0469, respectively. In Fig- 
ure Q] and Figure |2j we present G n)T for four different regimes when r takes values: 
0.020, 0.035, 0.045, 0.050, respectively. The values 0.020 and 0.035 correspond to 
'ultra' -sparse regime and emergence of a giant component, FigureQ] The values 0.045 
and 0.050 correspond to 'almost' -connected and connected regimes, Figure |2] 




Figure 1: 'Ultra' -sparse regime and the emergence of the giant component. 
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Figure 2: 'Almost' -connected and connected regimes. 

3 Bootstrap Percolation 

Bootstrap percolation (BP) is a cellular automaton defined on an underlying graph G = 
(V, E) with state space {0, 1} V whose initial configuration is chosen by a Bernoulli 
product measure. In other words, every node is in one of two different states or 1 
(inactive or active respectively), and a node becomes active with probability p inde- 
pendently of other nodes within the initial configuration. 

After drawing an initial configuration at time t = 0, a discrete time deterministic 
process updates the configuration according to a local rule: an inactive node becomes 
active at time t + 1 if the number of its active neighbors at t (not necessarily the nearest 
ones) is greater than some defined threshold 8. Once an inactive node becomes active 
it remains active forever. A configuration that does not change at the next time step is 
a stable configuration. A configuration is fully active if all its nodes of are active. 

An interesting phenomenon to study is metastability near a first-order phase transi- 
tion. Do there exist < p' c ^ p" c < 1 such that: 

Np < p') lim P„ (V becomes fully active) = , 

and 

(Vp > p") lim P p (V becomes fully active) = 1 ? 

Further, is it necessary for p' c to be asymptotically equal to p"? 

A study of BP on a regular infinite tree first appeared in j@J. Subsequently, the re- 
lations between the branching number of an infinite (non-regular) tree, threshold value, 
and p necessary to fully percolate the tree were studied in |5|]. 
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An example of BP is a d-dimensional lattice Z d equipped with Bernoulli product 
measure with 9 = d |Q]]. For Z d and V = [0, L — l] d the existence of a unique threshold 
p c was shown in Concretely for 1? and V = [0, L — l] 2 the exact threshold value 
is p c = 7r 2 /(181nL) ifnll . Furthermore the sharp threshold for bootstrap percolation 
in 7L d in all dimensions was provided in 0]. 

Additionally to BP on trees and lattices, there has been recent work of BP on ran- 
dom regular graphs Erdos-Renyi random graphs [19], as well as random graphs 
with a given degree sequence where the threshold depends upon node degree |2Q . 

3.1 Bootstrap Percolation on Connected RGGs 

The structure of G n ^ r is conducted by random positions of its nodes and radius r = 
r(n); so it is more 'irregular' than the structure of a tree or a lattice. In this work we 
are interested in BP on G n ^ r which for brevity we denote by BP{G n ^ r: p, 6). In this 
process a node becomes active with probability p independently of other nodes in the 
initial configuration and an inactive node becomes active at the following time step if 
at least 9 = 7D of its neighbors are active, where 7 = j(n) and D(n) = E(deg) = 
r 2 7r = a In n is the expected node degree. 

For the critical thresholds p' c and p" c in BP(G njr ,p, 0), we derive bounds p' ^ p' c 
and p" ^ p" c such that a connected G n r does not become fully active for p < p' whp, 
and conversely, becomes fully active for p > p" whp. These bounds are schematically 
presented in Figure [3] 

p' p " 7 

1 1 1 1 — 1 1 1 

p' c f (: 1 

Figure 3: Bounds (p',p") on the critical thresholds (p' c ,p"). 

The main ideas of the proofs are as follows. We obtain the distribution of the num- 
ber of active neighbors for each node at the initial configuration. For p < p' we use the 
Poisson tail bound and the union bound (see (fTOl ) in Appendix) to show that an initial 
configuration is stable whp. For p > p" we use the Bahadur-Rao theorem (see Claim|7] 
in Appendix) to lower bound the number of active neighbors for each node. Then 
we develop a geometric argument to show that a stable, fully active, configuration is 
reached within 0{y/n/r) steps whp. This geometric argument leverages the following 
simple observation about BP in 1? with = 1. 
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Figure 4: Functions H(x) and J(x). 

Lemma 1 Consider BP in Z 2 with the threshold 6 = 1 and the initial probability 
p > 0. For any N and p = uj(l/\fN), a square [0, iV] 2 becomes fully active within 
O(N) steps whp. 



We first introduce the following functions upon which our analysis will heavily 
depend. For the function H(x) := x In x — x + 1 on [0, +oo) (see Figure [4] left), define 
HjT 1 : [0,1] [0,1] to be the inverse of H(x) on [0,1], and : [0, +oo) -)• 

[1, +oo) to be the inverse of H(x) on [1, +oo). Analogously for the function J{x) := 
x^ 1 H(x) = Inx — 1 + x^ 1 on (0, +oo) (see Figure |4]right), define 1 : [0, +oo] — > 
[0, 1] to be the inverse of J(x) on [0, 1], and J^ 1 : [0, +oo) — > [1, +oo) to be the 
inverse of J(x) on [1, +oo). 

We now provide bounds on the critical thresholds on p in Theorem |2] and Theo- 
rem HJ 

Theorem 2 Consider bootstrap percolation BP(G ntr ,p,6) where r = V n~ 1 a In n 
and 6 = 7a In n. For a > 1, 7 € (0, 1) and when 

P<p' :=rf/J^ 1 (l/a'r), 
G ntr does not become fully active whp. 

Proof We show that for the conditions of the assertion, an initial configuration is 
stable. The number of active nodes in the initial configuration follows Poisson distri- 
bution Po(pn). The degree distribution of a node is Po(r 2 7r) — 1, and the expected 



degree D = r tt = a Inn. By the thinning theorem [22] the number of active neigh 



bors in the initial configuration follows Po(pD) — 1. Consider the activation rule in 
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BP(G njri p, 0). The probability that a node becomes active at the next time step given 
that it is inactive initially is P (Po(pD) — 1 ^ jD). 

For p > 7, given that »D — )■ oo, the tail bound on a Poisson random variable (fTOl ) 
implies for any node P (Po(pD) — 1 ^ 7-D) = 1 — o(l/n). Hence an initial configu- 
ration becomes fully active at the next time step with probability 1 — o(l). Therefore 
we consider the case p ^ 7 and seek a maximal p' ^ 7 (see Figure [3]> such that 
BP(G rhr ,p' , 7) does not become fully active whp. It follows 

P (Po(p£>) - 1 > 7D) < P (Po(pD) ^ jD) < exp(-pDH(j/p)) . (1) 

The same inequality (TTOb yields that the number of nodes Po(n) within the square 
[0, \fn] 2 is concentrated around its mean n whp. Hence the union bound over all nodes 
provides 

P p (the initial configuration is stable) ^ 1 — exp ((1 + o(l)) In n — pDH{^ /p)) . 

(2) 

Given D = a Inn, the condition paH{^/p) > 1 suffices that the initial configuration 
is stable whp. The function J(x) = x~ l H(x) is monotonically decreasing on (0, 1), 
monotonically increasing on (1, +00), with the minimum attained at x = 1. Hence 
for any positive 7 < +00 there are two solutions of J(x) = l/aj, denoted x\ < 
1 < X2- This yields p > 7/xi > 7 or p < 7/X2 < 7. The acceptable solution is 
p < 7/^2, since we consider the case p < 7. For J(^/p) > 1/07 from © it follows 
the probability that the initial configuration is stable tends to one as n tends to infinity. 
Finally, a bound on p is given by 

P<P' ^l/J^O-faj), 
which concludes the proof. □ 

The following result clarifies the feasible region for a and 7 in Theorem @] 
Lemma 3 The condition a ^ 5tt / H (5irj) is equivalent to: 

1 



7 G 



0,-^H^fa/a) 



, for a < to . 



and 
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^H^ 1 fa/ a), ^-H^ 1 fa /a) 



for a ^ to . 



Proof By inspection of the function H(x). □ 
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Theorem 4 Consider bootstrap percolation BP{G n>r ,p,6) where r = Vvr 1 a\nn, 
9 = ja Inn, and a > 1. Wzen a ^ 5ir/H(5ir^) and 7 G (0, l/5ir) for 
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p > p := mm < 7, 



57T7 



i ^(iM)j ' 

G ntr becomes fully active within 0{\/n/r) steps whp. 

For p > 7 an initial configuration becomes fully active at the next time step whp 
(see the proof of Theorem [2]), therefore we consider the case p ^ 7. 

The proof of Theorem[4]consists of two parts. We tile the square [0, y/n\ 2 into cells 
r/\/5xr/\/5 and show that in the initial configuration: (i)Whena > 1, 7 G (0,l/57r), 
and a ^ hit / H(hir^), every cell contains at least 7D nodes whp; (ii) When p > p" at 
least one cell contains 7Z? or more active nodes. By Lemma[T]it follows that for a, j,p 
in the specified ranges G n ^ r becomes fully active within 0(yjn/r) steps whp. 
Proof Tile the square [0, \/n\ 2 into cells r jxfh x r/v5, see Figure [5] Define the 
area of a cell A := r 2 /5 = alnn/57r. Call two cells neighboring if they share one 
side. Notice every pair of nodes within the same cell or within two neighboring cells 
are adjacent by the choice of the size of a cell. Define G' nr on the set of nodes of 
G n>r as follows. The set of edges of G' n r consists of the subset of edges of G n r whose 
terminal nodes belong to the same cell or two neighboring cells. Then the monotonicity 
of bootstrap percolation yields 

P p [G' nr becomes fully active) ^ P p (G n>r . becomes fully active) . (3) 

Therefore it is sufficient to show that whp G' n r becomes fully active when p > p" ([8]). 

Part (i) (To show that every cell contains at least 7D nodes whp.) We first bound 
the probability that an arbitrary cell contains at most jD nodes. The number of nodes 
in a cell follows Po(A), i.e., Po(alnre/57r). Moreover, the numbers of nodes in cells 
are independent random variables (given the Poisson point process X). For 7 ^ l/5ir, 
from (ITOb we obtain 

(/ Q In 71 \ 
Po ( — - — j 7a In n 

( a Inn 

s£ exp — H(5ttj) 

\ 07T 

= n 5ir v 11 , 

The total number of cells in [0, \/n} 2 is 5n/r 2 = 57rn/(olnn) = o(n). The union 
bound taken over all cells yields 

P (every cell contains at least -yD nodes) ^ 1 - o (n 1- ^ H(57r7) ) . (4) 
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Figure 5: Tiling the square [0, x/n]" 



Finally, for a ^ Bit / H , from (@]) it follows that every cell contains at least 
nodes whp. 

Part (ii). (To show that at least one cell contains 7D or more active nodes.) We now 
derive conditions such that at least one cell contains at least 9 = 7D active nodes in 
the initial configuration. In order to guarantee that whp there is at least one cell among 
5n/r 2 = O (n/ Inn), which contains at least 9 active nodes in the initial configuration, 
it suffices to find p such that 

(In ti \ 
— , (5) 

since 

lim 1 - (1 - w(lnn/n)) e(n/lnn) = 1 . 
Define a := 5n^/p — 1, then by rewriting ([5]) we need p such that 

p /ftKM)-M >n + J_\ M /W). (6) 

y pA pA J \ n J 

By the Bahadur-Rao tail bound (Hi, as n — > 00, i.e., pA = 0(plnn) — > 00, the shifted 
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Poisson random variable Po(pA) — pA satisfies 



S Po(pA)-pA ^ a | 1 
V pA pA 



vT+a 1 



exp (—pAI(a)) 



ay/27T y/pA 
where the rate function is defined by 

1(a) = sup {sa — e s + s + 1} = (1 + a) ln(l + a) — a . 
(See Appendix for the details.) Therefore for Q to be satisfied we require 



n 



VI + a 1 



Inn aV2K y/pA 
The left hand side of (O equals 



exp (-pAI(a)) = w(l) . 



(V) 



/fl «P r , ^, 3 1 pa VI + « 

exp 1 i (a) Inn m in n m h m ==- 

M 5vr V ') 2 2 5vr aV2T 



which is w(l) if 1 > apl(a)/bir. Given a = Bir^/p— 1, the condition 1 > apI(a)/hiT 
is equivalent to 1/07 > H(5n-f /p) / (bn^y /p), and moreover to 



p > p := 



57T7 



(8) 



To complete the proof notice that once any jD nodes within a cell become active, 
all nodes within that cell become active at the next time step as would all nodes within 
its neighboring cells. This resulting process which jointly activates all nodes within one 
cell is equivalent to activating a site in 1? . The resulting BP in 1? has the threshold 
9 = 1 by construction, see Figure [5] Thus BP(G' n r ,p, 9) becomes fully active when 
p > p" by LemmaQ] The proof follows from ([3]>. □ 



Remark 1 For non-trival percolation threshold, that is, p" ^ 7, it is necessary 

1 



0.55. 



Jr(5tt) 

Remark 2 When a > 1, the upper bound on 7 in Lemma\3\is further tightened: 

1 



7 € 



0. 



aJ R (57r) 



, for a < 57r , 



and 



7 e 



— tf^ 1 5vr/a , 

5f L aJ R (5n) 



, /or a ^ 57r . 
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3.2 Analysis of Bounds on Critical Thresholds 

The critical threshold p' = 7/ J^ 1 (l/cry) can be rewritten as 

lnp' = -lna-ln((l/a7)J fi 1 (l/a7)) . (9) 

The function — In (xJ^ l (x)) is monotonically decreasing in x, hence p' is mono- 
tonically increasing in a and monotonically decreasing in 7. As an example we numer- 
ically compute and tabulate p' for 7 = 1/20 and different values of a in Table 1. In 
Figure |6j p' is plotted as a function of a for different values of 

7 € {1/70, 1/60, 1/50, 1/40, 1/30, 1/20} . 



a 


p> 


p" 


3 


0.0000234198 


0.0003678767 


4 


0.0001242460 


0.0019516511 


5 


0.0003391906 


0.0053279940 


6 


0.0006649716 


0.0104453500 


7 


0.0010794693 


0.0169562642 


8 


0.0015576467 


0.0244674579 


9 


0.0020779022 


0.0326396121 


10 


0.0026234549 


0.0412091329 


25 


0.0101188498 


0.1589465210 


50 


0.0174952121 


0.0174952120 


100 


0.0246619916 


0.3873896589 



Table 1: Bounds p',p" on the critical thresholds for different values of a when 7 = 
1/20. 

The experiments are performed on G n ,r with n = 15000 and n = 25000 nodes, 
and r = ^/alnn/iT for the cases: (i) a = 30 and 7 = 1/100, and (ii) a = 35 and 
7 = 1/75. On these instances of graphs, for each chosen value of p in (0, 1) we 
simulate BP 100 times. More precisely, within each experiment we generate a random 
initial configuration with the probability p and perform BP with the threshold 9 = jD 
where the expected degree D is calculated for a given input G n ^ r . 

Numerical results are presented with the initial probability p on the horizontal axis, 
and the percentage of fully active stable configurations on the vertical axis. Four cases 
when (a = 30, 7 = 1/100), (a = 35, 7 = 1/75), for n = 15000, 25000, are presented 
in Figures 13 [U|9] and [lOl respectively. These charts match the bounds derived theoret- 
ically for p' and p". Further, they appear to support the case that p' c ^ p" c even though 
we do not currently have a proof one way or the other. 
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Figure 7: Percentage of fully percolated configurations in 100 simulations of 

BP(G n . r ,p,6) when a = 30,7 = 1/100, n = 15000, r = ^30 lnn/vrn rj 0.07824, 
D = 30 Inn rj 288.47 and 6» = [100- 1 E(deg)] = [2.88] = 3. The bounds are 
p' = 0.000133 and p" = 0.002089. 
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Figure 8: Percentage of fully percolated configurations in 100 simulations of 

BP(G n ,r,p,6) when a = 30,7 = 1/100, n = 25000, r = ^30 lnn/vrn » 0.06219, 
D = 30 Inn « 303.80 and 6> = [100" 1 E(deg)] = [3.04] = 4. The bounds are 
p* = 0.000133 and p" = 0.002089. 
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Appendix 

Lemma 5 ( Concentration on a Poisson random variable, see KlAl ). A Poisson random 
variable Po(A) (with X > 0) satisfies: 

P(Po(A)^A;) < exp (—XH(k/X)) , for k^ X, (10) 
P(Po(A)<fc) ^ exp(-XH(k/X)) , for k < A, (11) 

where H(x) = x In x — x + I for x > 0. 



Theorem 6 (Bahadur-Rao, see ATj/J Let X±,X2, ... be an Ltd. sequence of random 
variables such that E(ATi) = and M(s) := K(e sXi ) < oo for all sei If X\ is of 
lattice type and ¥(Xi = a) > 0, then 

lim P ( 4r V Xi ^ a] exp (NI(a)) \fN = — = , 

\ N h. ) ay/to (l-exp(-8 a ))' 

where 1(a) := sup sgR (sa — lnM(s)) attained at s = s a , and a 2 = M"(s a )/M(s a ) — 



a 2 . 



Claim 7 A Poisson random variable Po(N)for N — > oo satisfies 



lim p ( Po(N) ~ N A exp ^ = VT+a 

where 1(a) = (1 + a) ln(l + a) — a for a ^ 0. 
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Proof Let X{ ~ Po(l) — 1 for i = 1,2,... be the independent lattice type ran- 
dom variables. We have E(JQ) = and Var(Xj) = 1. Consider: (i) the mo- 
ment generating function M(s) := K(e sXi ) = exp(e s — s — 1), (ii) the rate func- 
tion 1(a) = sup sGR (sa — lnM(z)) = (1 + a) m(l + a) — a which is attained at 
s a = ln(l + a), and (iii) the variance a 2 := M"(s a )/M(s Q ) — a 2 = 1 + a. Now the 
claim follows from Theorem [6] □ 
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